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Abstract 

We have derived constraints on the tau neutrino mass and fourth generation 
mixing from an analysis of the partial widths of tau lepton decays, in particular: 
T~^e~Ue^T, — >yU~z^^z^T, T—^'7i~h'r, and T—^K'^r- We present predictions 
for the tau decay widths, allowing for a non-zero tau neutrino mass, rriiy^ , and 
for mixing with a neutrino of mass m^^ > Mz/2, which is parametrised using a 
Cabibbo-like mixing angle, 6'l. By comparison of these theoretical predictions 
with the experimental measurements, we obtain the following bounds at the 
90% confidence level: m^,^ < 42 Me F and sin^ Oi, < 0.014. 



1 Introduction 



Massive neutrinos feature in many extensions to the Standard Model They have been suggested 
as a source of dark matter and, if they oscillate in the sun, as an explanation for the deficit in the 
observed solar neutrino flux [^Q. The best experimental upper limit on the tau neutrino mass is 
nii,^ < 24 MeV at the 95% confidence level , which was obtained using many-body hadronic decays 
of the r. 

In this paper we describe a complementary method for constraining the r neutrino mass, m,^^, 
from precise measurements of r partial widths for the following decays[| :r~ — i>e~i^eZ/^, t~ ^ 

ti'i't, and r~ K'v^. The dependence of the purely leptonic decay rates on rriy^ has been 
considered by others f^-H whereas, to the best of our knowledge, the hadronic decays have not 
previously been analysed for this purpose. The r partial decay widths are also sensitive to mixing 
in the leptonic sector with a fourth generation weak isospin lepton doublet (i/l,L~) P,|9|-[13|, where 
the neutrino has a mass M > Mz/2. 

We present calculations of the r partial widths for these channels which allow for the effects of non- 
zero tau neutrino mass and mixing. Then, using data from e^e~ experiments, we derive constraints 
on the tau neutrino mass and mixing. 

This approach is complementary to traditional analyses of the kinematic end-point of multi-hadron 
tau decays, which yield tighter constraints on m,^^ but which are insensitive to fourth generation 
mixing. Moreover, the channels we consider are statistically independent and theoretically very well 
understood. 



2 Partial widths for r decays 



In the Standard Model describing the electroweak interaction, the partial width Ti for the decay 
r"— >£"z/^z/^(Xem), with i~ = e~ , fi~ and Xem = 1, 11, e+e", . . ., is given by 



where Gp is the Fermi constant, vtir and are the r mass and lifetime. The radiative-correction 



function Rg has been calculated [p^-[17|, with ainir) — 1/133.3 |T^, to be 
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If the neutrino masses are zero for all generations, then the phase-space factor, F^, is given by the 
well-known expression: 



Fj^(x) = l-8x- UxHnx + 8x 



1, (£ = e); 

0.9726, {i = iJ,); 



(3) 



where x = ml/ml. A value of Gp = (1.16639 ± 0.00002) x lO^^GeV^^ ^ jg obtained from the 
measurement of the muon lifetime using equation ||, with £^ = e^ and substituting r — ^ /i. Gp 
implicitly includes the residual effects of radiative corrections not explicitly included in equation 0. 



-'^ Henceforth we denote the branching ratios for these processes as Be, B-j^^ Bk respectively; Bi denotes either Be 
or Bfj, while Bh denotes either B^^ or Bk- 



The partial widths for the decays r ^h i^ri'l), with h = it /K , are given by 

Th = ^ = ^flR,Kp\'F^, (4) 

where /h are the hadronic form factors, /^r and /k, and V^^ are the CKM matrix elements, V^d 
and Vus, for vr" and K~ respectively. From an analysis of vr" — /U~z/^ and K~ fJ'~'^n decays, one 
obtains A|Kd| = (127.4 ± 0.1)MeV and /k|Ks| = (35.18 ± 0.05)MeV P, and references therein]. 
The radiative correction factor, i?h, is given by 

i?, = l + ^lnf!^)+...^1.02. (5) 

The ellipsis represents terms, estimated to be (9(±0.01) |l^, which are neither explicitly treated nor 
imphcitly absorbed into Gp, /7r|Kid|, or /k|Kis|- For massless neutrinos the phase-space function, 
is given by 

^o.x^^_x2^f 0.9877, (h = 7r); 
^hW U X) - 1 0.8516, (h = K); 

where mh is the hadron mass and x = m\/m1. 

We now consider the effects of neutrino masses and mixing on the r decay rates. The expression 
for r^, allowing for non-zero neutrino masses and mixing between n lepton generations, is given by: 

^^ = ^,R^tt\Urm^'F,, (7) 
-L^^" i=l j=l 



where Uai is the lepton mixing matrix |]T8|, page 276] and F^, the phase-space factor, depends on the 
neutrino and charged lepton masses. 



The electron and muon neutrinos couple dominantly to vi and V2 and have small masses |T8|. We 
therefore neglect the masses of vi and V2 and their mixing with v^. We do, however, consider the 
possible existence of a fourth generation neutrino, i^l, of mass niy^ > Mz/2 which mixes with z/,- (LEP 
data have constrained the number of neutrinos of mass m^, < Mz/2 to be N^y = 2.991 ± 0.016 ||18|| ). 
Since such a neutrino is kinematically forbidden in r decays, the corresponding phase-space factor is 
zero, is therefore suppressed by a factor which depends on the strength of the mixing of the third 
and fourth generations, but not on the mass m^^. From the above, equation |^ reduces to 

r, = j^Ml- sin' 9^)F„ (8) 

where the factor of (1 — sIu^^^l) = |f/r3p allows for the Cabibbo-like suppression due to fourth 
generation mixing, and the phase-space factor is given by: 

F,{x, y) = F^{x) - %y{l - xf + ■ ■ ■ (9) 

where y = ml^/m'^, Ff is given by equation |^ and the ellipsis represents negligible terms of higher 
order in iriy^. Figures 1(a) and 1(b) show the variation of FijF^ with vtiy,^ for — ^^v^v-y and 
r~ — s>/i~z/^z/^ decays respectively. 

The expression for Fh, allowing for non-zero neutrino masses and mixing in a similar fashion, is 
given by: 

Th = fhRh\Vag\'^y^AUri\ F^ 

-"-"^ i=l 

= ^flR^\Vapni-sm'e^)F^ (10) 



where the phase-space function has been calculated to be: 



F^{x,y)=F^ 
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with y = ml.^/m'^. Figures 1(c) and 1(d) show the variation of F^/F^ with rriiy^ for T^n^Ur and 
r— >K~z/t- decays respectively. Note the lower sensitivity to mj/3 of the (two-body) hadronic decay 
modes compared to the (three-body) leptonic modes, despite the higher masses of the final state 
hadrons. 



3 Constraints on mass and mixing 

In the standard model of electroweak interactions the three lepton generations interact in an identical 
way, apart from effects caused by their differing masses. In particular, they each couple with the 
same strength to the charged weak current. If this assumption of universality is relaxed for the 
r, then the G| factor appearing in equations |1| and ^ may be replaced by G^Gp, where Gp is not 
necessarily equal to Gp. If the measured values of Gp, evaluated assuming massless neutrinos and 



no mixing, appeared to be significantly smaller than Gp it could indicate either new physics [13] or 
the suppression of decay rates due to non-zero neutrino masses or mixing. 



We use the Particle Data group values and errors |T^ for the measured quantities, in particular: 
Tr = (291.0 ± 1.5)fs; rrir = (1776.96l|ji?t|j:?^)MeV; 

Be = (17.80 ±0.08)%; = (17.30 ± 0.10)%; 

B^ = (11.07 ±0.18)%; Bk = (0.71 ± 0.05)%. 

The sole exception is the use of the BES r mass obtained from the measurement of r+r^ 
production near threshold, since this has no dependence on m^^ . Substituting in equations |1| and ^ 



for the measured quantities, and assuming that m^^ = and sin = we obtain 

r 0.996 ±0.007 
G£ _ I 0.996 ±0.008 
Gv ~ I 1.010 ±0.017 
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where the number of standard deviations from the universality prediction (of unity) are shown in 
parentheses. These results, which are all consistent with unity, indicate that the lepton couplings 
are universal and show no indications for non-zero neutrino mass or mixing. We therefore assume 
universality holds and use the measured r partial widths to constrain m^,^ and sin^^L. 

We have mapped the likelihood of observing the measured r partial widths, as a function of both 
rriu^ and sin^ ^l, by randomly sampling all the quantities used according to their experimental errors, 
allowing for the estimated 1% theoretical uncertainty on R^. The CLEO measurement of the r mass 
was used to further constrain rriu.^. From an analysis of r+r~ — (vr'^nTr^z/^) (7r~m7r°z/^) events (with 
n < 2,m < 2,1 < n + m < 3), CLEO determined the r mass to be rrir = (1777.8 ± 0.7 ± 1.7) ± 
[m^^{MeV)]'^/UOOMeV 0. The likehhood for the CLEO and BES measurements to agree, as a 
function of m^^ is included in the global likelihood. 

Figure 2(a) shows the 90% and 95% contours of the two dimensional likelihood distribution com- 
bined for all four r decay channels. No evidence is seen for a non-zero neutrino mass, nor for 
mixing. By integration of the two dimensional likelihood over all values of sin^ we obtain the 
one-dimensional likelihood for m^^, independent of the value of sin^^L, as shown by the solid line 



of figure 2(b). We obtain the following upper limits: 771^,3 < 42(48) Me Fat the 90(95)% confidence 
levels. The solid line of figure 2(c) shows the one-dimensional likelihood distribution for sin^^^L, in- 
tegrated over all values of rriu.^, from which we derive the upper limits: sin^ < 0.014(0.017) at the 
90(95)% confidence levels. Table 1 summarises the upper limits of ruy^ and sin^6'L, individually for 
each channel and for all channels combined. Since the mixing of 777.^,3 with other neutrinos is small, 
the limits derived for 777^3 can be reasonably interpreted as limits on 777,^^. 

Non-inclusion of the constraint from the CLEO r mass measurement results in the following limits: 
771^3 < 59(68) Mel/ and sin^ < 0.013(0.016) at the 90(95)% confidence levels. The slight reduction 
in the sin^6'L limit values, despite the increase in the 772.1^3 limit, is due to the anticorrelation of the 
effects on the hkelihood of non-zero values for 777,^3 and sin^ ^l- 

4 Conclusions 

We have derived the following constraints on the tau neutrino mass and fourth generation mixing from 
an analysis of the partial widths of tau lepton decays: 777,^^ < 42 Me V and sin^ < 0.014 at the 90% 
confidence level. These resiilts are statistically and systematically independent of traditional end- 
point analyses using multi-liadronic decays of the r. Moreover we simultaneously consider mixing, 
and are insensitive to fortuitous or pathological events close to the kinematic limits, details of the 
resonant structure of multi-hadron r decays, and the absolute energy scale of the detectors. 

These results will improve with more precise measurements of the r mass, lifetime, and branching 
fractions, for example at a r-charm factory. Ultimately, we expect that this method will be limited 
by the uncertainty on the r lifetime. The extension of this technique to include multi-hadronic r 
decays, in conjunction with an improved theoretical description, should provide considerably more 
sensitivity due to the higher multiplicity of the final states. 
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Table captions 



Table 1. 

Upper limits obtained on 171,^^ and sin^ 6'l, seperately for each r decay channel and for all channels 
combined. 



Figure captions 

Figure 1. 

Relative phase-space suppression factors, as a function of tau neutrino mass for (a) t~ -^e~Ue^T, (b) 
T~ ^ /i'PhI't-, (c) t~— >7r~i/^, and (d) r~— >K"i/^. 



Figure 2. 

Likelihood distributions for all r decay channels combined, for (a) 171,^^ vs. sin^ 6'l, (b) 171,^^, integrated 
over sin^^L, and (c) sin^^L, integrated over 171^3 ■ 
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Abstract 

We have derived constraints on the tau neutrino mass and fourth generation 
mixing from an analysis of the partial widths of tau lepton decays, in particular: 
T~ — )• e~UeV-r, T~ — )• fi~u^u-r, T — )• Tr~u-r, and T — )• K~u-r. We present predictions 
for the tau decay widths, allowing for a non-zero tau neutrino mass, m,^^, and 
for mixing with a neutrino of mass m,^ > Mz/2, which is parametrised using a 
Cabibbo-like mixing angle, 6i,. By comparison of these theoretical predictions 
with the experimental measurements, we obtain the following bounds at the 
90% confidence level: m^^ < 42 Me and sin^ < 0.014. 



1 Introduction 



Massive neutrinos feature in many extensions to the Standard Model [1]. They have been suggested 
as a source of dark matter [2] and, if they oscillate in the sun, as an explanation for the deficit in the 
observed solar neutrino flux [3,4]. The best experimental upper limit on the tau neutrino mass is 
m,y^ < 24MeV at the 95% confidence level [5], which was obtained using many-body hadronic decays 
of the T. 

In this paper we describe a complementary method for constraining the r neutrino mass, m,^^, 
from precise measurements of r partial widths for the following decays^ :r~ — )• e~z/ei'T, t~ ^ fi~u^u-r, 
T~ — )• Tr~u-r, and r~ — )• K~u-r. The dependence of the purely leptonic decay rates on m,^^ has been 
considered by others [6-8] whereas, to the best of our knowledge, the hadronic decays have not 
previously been analysed for this purpose. The r partial decay widths are also sensitive to mixing 
in the leptonic sector with a fourth generation weak isospin lepton doublet (z/l,L~) [6,9-13], where 
the neutrino has a mass M > Mz/2. 

We present calculations of the r partial widths for these channels which allow for the effects of non- 
zero tau neutrino mass and mixing. Then, using data from e"'"e~ experiments, we derive constraints 
on the tau neutrino mass and mixing. 

This approach is complementary to traditional analyses of the kinematic end-point of multi-hadron 
tau decays, which yield tighter constraints on m,^^ but which are insensitive to fourth generation 
mixing. Moreover, the channels we consider are statistically independent and theoretically very well 
understood. 



2 Partial widths for r decays 

In the Standard Model describing the electroweak interaction, the partial width for the decay 
T~ £~ i'iUr{XEM) 1 with £~ = e~,/i~ and Xem = 7, 77, e+e", . . ., is given by 

where Gy is the Fermi constant, and are the r mass and lifetime. The radiative-correction 
function Ri has been calculated [14-17], with a(mT-) ~ 1/133.3 [17], to be 



27r V 4 / 



3 m 



2 



5 



0.9960. (2) 



If the neutrino masses are zero for all generations, then the phase-space factor, is given by the 
well-known expression: 

= 1 - 8. - nA. + 8.= - c { : (3) 

where x = ml/ml. A value of Gp = (1.16639 ± 0.00002) x lO-^GeV-^ [18] is obtained from the 
measurement of the muon lifetime using equation 1, with i~ = e~ and substituting t ^ fi. Gp 
implicitly includes the residual effects of radiative corrections not explicitly included in equation 2. 



^Henceforth we denote the branching ratios for these processes as B^, Bf^, 3^^, Bj^ respectively; Bi denotes either B^ 
or while Bh denotes either or Bk- 



1 



The partial widths for the decays r — )-h i't(7), with h = tt /K , are given by 

T^^^ = ^flR^MF^, (4) 

where /h are the hadronic form factors, and /k, and are the CKM matrix elements, V-ad 
and V-as, for 7r~ and K~ respectively. From an analysis of 7r~ — )• /i~z/^ and K~ — )• /i~z/^ decays, one 
obtains /^|Kd| = (127.4 ± 0.1)MeV and /k|Ks| = (35.18 ± 0.05)MeV [19, and references therein]. 
The radiative correction factor, i2h, is given by [17] 

i2H = l + -lnf^)+...^1.02. (5) 

The ellipsis represents terms, estimated to be (9(ib0.01) [19], which are neither explicitly treated nor 
implicitly absorbed into Gp, f-wlVadl^ or /k|Kis|- For massless neutrinos the phase-space function, F°, 
is given by 

J^i,{x)-[L X) - |Q_85;Lg^ (h = K); 

where mn is the hadron mass and x = m^/m^. 

We now consider the effects of neutrino masses and mixing on the r decay rates. The expression 
for r^, allowing for non-zero neutrino masses and mixing between n lepton generations, is given by: 

1=1 j=i 

where Uai is the lepton mixing matrix [18, page 276] and F^, the phase-space factor, depends on the 
neutrino and charged lepton masses. 

The electron and muon neutrinos couple dominantly to Ui and U2 and have small masses [18]. We 
therefore neglect the masses of Ui and U2 and their mixing with U3. We do, however, consider the 
possible existence of a fourth generation neutrino, z^l, of mass m,^ > Mz/2 which mixes with u-r (LEP 
data have constrained the number of neutrinos of mass m,^ < Mz/2 to be N,^ = 2.991 it 0.016 [18]). 
Since such a neutrino is kinematically forbidden in r decays, the corresponding phase-space factor is 
zero, is therefore suppressed by a factor which depends on the strength of the mixing of the third 
and fourth generations, but not on the mass m,^. From the above, equation 7 reduces to 

r. = ^i2.(l-sin^W, (8) 

where the factor of (1 — sin^ ^l) = iC^rSp allows for the Cabibbo-like suppression due to fourth 
generation mixing, and the phase-space factor is given by: 

F^{x,y) = F^{x)-8y{l-xY + --- (9) 

where y = ml^/m^, F° is given by equation 3 and the ellipsis represents negligible terms of higher 
order in m,^^. Figures 1(a) and 1(b) show the variation of Fi/F° with m,^^ for r~ — )• e~UeV-r and 
T~ ^ fi~u^u-r decays respectively. 

The expression for Fh, allowing for non-zero neutrino masses and mixing in a similar fashion, is 
given by: 

Glm^ 



%^/hi2h|K/3r(l-sin^W (10) 



where the phase-space function has been calculated to be: 



F^{x,y) = F^{x) 



2 + x-y 



\ 



2 + 2x-y 



(1 



xy 



(11) 



with y = ml^/m^. Figures 1(c) and 1(d) show the variation of Fy^j with 7X1^,^ for T^Tr~u-r and 
r— )-K~z/t- decays respectively. Note the lower sensitivity to m,^^ of the (two-body) hadronic decay 
modes compared to the (three-body) leptonic modes, despite the higher masses of the final state 
hadrons. 



3 Constraints on mass and mixing 

In the standard model of electroweak interactions the three lepton generations interact in an identical 
way, apart from effects caused by their differing masses. In particular, they each couple with the 
same strength to the charged weak current. If this assumption of universality is relaxed for the 
r, then the Gp factor appearing in equations 1 and 4 may be replaced by GpGp, where Gp is not 
necessarily equal to Gp- If the measured values of Gp, evaluated assuming massless neutrinos and 
no mixing, appeared to be significantly smaller than Gp it could indicate either new physics [13] or 
the suppression of decay rates due to non-zero neutrino masses or mixing. 

We use the Particle Data group values and errors [18] for the measured quantities, in particular: 

= (291.0 ± 1.5)fs; = (1776.96iSi?iS:?^)MeV; 

Be = (17.80 ±0.08)%; = (17.30 ± 0.10)%; 

= (11.07 ±0.18)%; Bk = (0.71 ± 0.05)%. 

The sole exception is the use of the BES r mass [20], obtained from the measurement of r"'"r~ 
production near threshold, since this has no dependence on m,^^. Substituting in equations 1 and 4 
for the measured quantities, and assuming that m^^ = and sin^ = we obtain 



Gp 



0.996 ± 0.007 (-0.54(7), for r" 
0.996 ± 0.008 (-0.55(7), for r" 



Gf ] 1.010 ±0.017 (±0.60(7), for r-^TT-z/^; '^^^'^ 
[ 0.987 ± 0.070 (-0.18(7), forr-^R-z/^. 

where the number of standard deviations from the universality prediction (of unity) are shown in 
parentheses. These results, which are all consistent with unity, indicate that the lepton couplings 
are universal and show no indications for non-zero neutrino mass or mixing. We therefore assume 
universality holds and use the measured r partial widths to constrain m^^ and sin^ 0i,. 

We have mapped the likelihood of observing the measured r partial widths, as a function of both 
m^^ and sin^ 0i,, by randomly sampling all the quantities used according to their experimental errors, 
allowing for the estimated 1% theoretical uncertainty on i2h- The CLEO measurement of the r mass 
was used to further constrain m,^^. From an analysis of r"'"r~ — )• (ir'^mr^u-r) {ir~rmr°u-r) events (with 
n < 2,m < 2, 1 < n ± m < 3), CLEO determined the r mass to be = (1777.8 ± 0.7 ± 1.7) ± 
['m^,{MeV)Y/UOOMeV [21]. The likelihood for the CLEO and BES measurements to agree, as a 
function of m,^^ is included in the global likelihood. 

Figure 2(a) shows the 90% and 95% contours of the two dimensional likelihood distribution com- 
bined for all four r decay channels. No evidence is seen for a non-zero neutrino mass, nor for 
mixing. By integration of the two dimensional likelihood over all values of sin^ 0i, we obtain the 
one- dimensional likelihood for m^^, independent of the value of sin^^L, as shown by the solid line 
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of figure 2(b). We obtain the following upper limits: m^^j < 42(48) Me at the 90(95)% confidence 
levels. The solid line of figure 2(c) shows the one- dimensional likelihood distribution for sin^ 0i,, in- 
tegrated over all values of m^^, from which we derive the upper limits: sin^ 0i, < 0.014(0.017) at the 
90(95)% confidence levels. Table 1 summarises the upper limits of m^^ and sin^^L, individually for 
each channel and for all channels combined. Since the mixing of m,^^ with other neutrinos is small, 
the limits derived for m,^^ can be reasonably interpreted as limits on m,^^. 

Non-inclusion of the constraint from the CLEO r mass measurement results in the following limits: 
m^3 < 59(68) Me and sin^ < 0.013(0.016) at the 90(95)% confidence levels. The slight reduction 
in the sin^ 0i, limit values, despite the increase in the m^^ limit, is due to the anticorrelation of the 
effects on the likelihood of non-zero values for m^^ and sin^ 0i,. 

4 Conclusions 

We have derived the following constraints on the tau neutrino mass and fourth generation mixing from 
an analysis of the partial widths of tau lepton decays: m^^ < 42 Me and sin^ 0i, < 0.014 at the 90% 
confidence level. These results are statistically and systematically independent of traditional end- 
point analyses using multi-hadronic decays of the r. Moreover we simultaneously consider mixing, 
and are insensitive to fortuitous or pathological events close to the kinematic limits, details of the 
resonant structure of multi-hadron r decays, and the absolute energy scale of the detectors. 

These results will improve with more precise measurements of the r mass, lifetime, and branching 
fractions, for example at a r-charm factory. Ultimately, we expect that this method will be limited 
by the uncertainty on the r lifetime. The extension of this technique to include multi-hadronic r 
decays, in conjunction with an improved theoretical description, should provide considerably more 
sensitivity due to the higher multiplicity of the final states. 
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Table captions 



Table 1. 

Upper limits obtained on m^^ and sin^^L, seperately for each r decay channel and for all channels 
combined. 



Figure captions 

Figure 1. 

Relative phase-space suppression factors, as a function of tau neutrino mass for (a) r~ -^e~UeV-r, (b) 
T~ ^ fi~u^u-r, (c) r~ — )• 7r~z/T-, and (d) r~— )-K~z/t-. 



Figure 2. 

Likelihood distributions for all r decay channels combined, for (a) m^^ vs. sin^ 0i,, (b) m^^, integrated 
over sin^ 0i,, and (c) sin^ 0i,, integrated over m^^. 
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Upper limit on m,^^ (MeV) 
at the 90% (95%) C.L. 


Upper limit on sin^ 9i, 
at the 90% (95%) C.L. 


T — )• e UeU-r 
T~ — )• Tr~U-r 


63 (72) 
66 (76) 
141 (160) 
333 (375) 


0.016 (0.019) 
0.018 (0.022) 
0.028 (0.033) 
0.036 (0.038) 


All channels 


42 (48) 


0.014 (0.017) 



Table 1. 
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